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A stellar prominence in the white dwarf/red dwarf binary 
QS Vir: evidence for a detached system 
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ABSTRACT 

Using high resolution UVES spectra of the eclipsing Post Common Envelope Binary 
QS Vir we detect material along the line of sight to the white dwarf at orbital phase 
4> = 0.16. We ascribe this to a stellar prominence originating from the M dwarf sec- 
ondary star which passes in front of the white dwarf at this phase. This creates sharp 
absorption features in the hydrogen Balmer series and Ca II H and K lines. The small 
size of the white dwarf allows us to place tight constraints on the column density of 
hydrogen in the n = 2 level of log 10 N2 — 14.10 ± 0.03 cm -2 and, assuming local ther- 
modynamical equilibrium, the temperature of the prominence material of ~ 9000K. 
The prominence material is at least 1.5 stellar radii from the surface of the M dwarf. 
The location of the prominence is consistent with emission features previously inter- 
preted as evidence for Roche lobe overflow in the system. We also detect Mg 11 448lA 
absorption from the white dwarf. The width of the Mg 11 line indicates that the white 
dwarf is not rapidly rotating, in contrast to previous work, hence our data indicate 
that QS Vir is a pre-cataclysmic binary, yet to initiate mass transfer, rather than a 
hibernating cataclysmic variable as has been suggested. 
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1 INTRODUCTION 

A binary system initially composed of two main-sequence 
stars, in an orbital period of between ~ 10- 1000 days 
jKolb fc Stehklll996h . will result in the two stars orbiting 
within a common envelope of material when the more mas- 
sive primary star evolves off the main-sequence and expands 
beyond its Roche lobe. The effect of this common envelope 
phase is to reduce the orbital period of the system dramati- 
cally. The resulting system, now composed of a tightly bound 
white dwarf and low-mass main-sequence star, evolves to- 
wards a semi-detached state via the loss of orbital angular 
momentum. Eventually, the surface of the main-sequence 
star will touch its Roche lobe and initiate mass transfer on 
to the white dwarf via a gas stream that passes through the 
LI point, creating a cataclysmic variable (CV) system. 

The evolution before and during the CV stage is 
drive n by angular momentum loss v ia gravitational radia- 
tion (IKraft et al . 1962; Faulk nerll971h and magne tic braking 
jVerbunt fc Zwaanlll98ll ; iRappaport et alJll983h . However, 
whilst the angular momentum loss from gravitational radi- 
ation is well understood, angular momentum loss rates via 
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magnetic braking are still uncertain and only recently has 
it become possible potentially to measure the angular mo- 
mentum loss r ate directly from pos t common envelope bina - 
ries (PCEBs) jParsons et alj|2010h . ISchreiber et atl l|201Ch 
recently found evidence of disrupted magnetic braking in 
PCEBs. 

Mass loss from nova eruptions during the CV stage 
are expected to result in a widening of the binary sepa- 
ration and an increase in the Roche radius of the donor 
star l|Shara et alj|l986h . This can result in the donor losing 
contact with its Roche lobe and thus mass transfer ceases 
and the systems becomes detached once more, until angular 
momentum once more brings the two stars close enough to 
re-initiate mass transfer. The existence of these hibernat- 
ing CVs is controversial, however, if they exist they would 
resemble detached pre-cataclysmic binaries. 

QS Vir (EC 13471 - 1258) is a close eclipsing binary 
system containing a DA white dwarf and a dMe dwarf. 
lO'Donoghue et al. (2003) presented photometric monitoring 
and time-resolved optical and ultra-violet spectroscopy of 
QS Vir and proposed that the system was a hibernating CV 
rather than a pre-CV. Their evidence for this was the large 
rotational velocity of the white dwarf (V ro t,WD = 400 ± 100 
kms -1 ), which implied that there had been an earlier phase 
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Table 1. Journal of VLT/UVES spectroscopic observations. 



Date Start Exp Orbital Conditions 

(UT) time (s) phase (Transparency, seeing) 

2001/08/19 00:07 300.0 0.16 Good, - 1 arcsec 

2002/04/23 01:45 300.0 0.05 Good, - 1.5 arcsec 

2007/08/17 00:46 1800.0 0.41 Variable, ~ 2 arcsec 



3 RESULTS 

3.1 Spectral features 

Figure [T] shows the 2001 spectra. At short wavelengths, 
where the white dwarf dominates, broad Balmer absorp- 
tion from the white dwarf is present. Additionally, hydrogen 
Balmer, Ca II H and K and Si I 3906A emission from the red 
star is evident. The Si I emission is irradiation driven there- 
fore it appears quite narrow in the 2001 spectrum since the 
heated face of the secondary star is not fully visible at the 
orbital phase that this spectrum was taken. At long wave- 
lengths absorption features from the red star dominate the 
spectrum. Similar features are seen in the 2002 and 2007 
spectra. 

The 2001 spectra also contain sharp absorption features 
corresponding to the hydrogen Balmer series and Ca n H and 
K lines; these are labelled in Figure [1] These sharp absorp- 
tion features are not present in the 2002 or 2007 spectra. We 
also detect Mg n 448lA absorption in all the spectra (see 
section [5^21 for a discussion of this feature). 

3.2 The origin of the sharp absorption features 

The sharp absorption features, labelled in FigurefT] are seen 
only in the 2001 spectra which means that they are either 
a transient feature or only visible at specific orbital phases. 
The depth of these features indicate that they cannot origi- 
nate from the photosphere of the red star because they are 
deeper than the total contribution from the red star, there- 
fore it is light originating from the white dwarf which is be- 
ing absorbed. Hence this absorption is from material cooler 
than the white dwarf, between us and the white dwarf at an 
orbital phase of 0.16. 

If this absorption is caused by a disk of material orbiting 
the white dwarf then the small change in direction of the 
material as it crosses the face of the white dwarf (slightly 
blueshifted on one side and slightly redshifted on the other) 
will broaden the absorption lines. Hence the width of the 
narrowest line, AV, can be used to estimate the minimum 
distance of the material from the white dwarf, R, we find 

where K scc is the projected orbital velocity of the M dwarf 
(266 kms^ 1 ), q = M sec /M-wt>, the mass ratio (0.52), a is 
the orbital separation (1.28-Rq) and -Bw d ~ 0.01 1J?0 is the 
radiu s of the white dwarf (all values from lO'Donoghue et al.l 
2003). The width of the narrowest line is 18 kms 1 (see 
Section 14. ip hence the distance of the material from the 
white dwarf is at least R > .22a. Given the relativel y low 
inclination of QS Vir (~ 75° , lO'Donoghue et al.ll2003T l it is 
unlikely that these absorption features are the result of disk 
material since it would have to be located far from the white 
dwarf and far out of the orbital plane. 

The velocities of the absorption features are consistent 
with the white dwarf's radial velocity (see Section ^. ll for ra- 
dial velocity measurements). This implies that the material 
is locked in rotation with the binary since any material ro- 
tating within the system, along the line of sight of the white 
dwarf, will have the same radial velocity as the white dwarf. 

based on these arguments the most likely origin of this 



of spin-up by mass transfer, and the behaviour of the multi- 
component Ha emission-line profile, which they concluded 
showed evide n ce of a weak mass transfer steam. Recently 
iRibeiro et all (|2010h questioned this conclusion, arguing 
that the emission feat ure seen close to the whit e dwarf in 
the Doppler images of lO'Donoghue et al.l ([2003) is the re- 
sult of wind accretion onto the white dwarf rather than a 
result of a low-level accretion stream shock. Their data also 
provided no indication that QS Vir was a hibernating CV. 
It is important to establish the evolutionary status of QS 
Vir since its eclipsing nature makes it ideal to study angular 
momentum loss. 

Here we present high-resolution UVES spectra of QS 
Vir and use them to identify and characterise a stellar promi- 
nence. We then discuss the role of the prominence within the 
system and evolutionary status of QS Vir as a whole. 



2 OBSERVATIONS AND THEIR REDUCTION 

Three spectra of QS Vir were obtained using the Ultravi- 
olet and Visual Echelle Spectrograph (UVES) installed at 
the European Southern Observatory's Very Large Telescope 
(ESQ VLT) 8 . 2-m t elescope unit on Cerro Paranal in Chile 
i Dekker et afll200Ch within both the 167.D-0407 (2001 and 
2002) and 079.D-0276 (2007) programs, details of these ob- 
servations are listed in Table Q] 

The reduction of the raw frames was conducted using 
the most recent standard pipeline release of the UVES Com- 
mon Pipeline Library (CPL) recipes (version 5.0.1) within 
ESORex, the ESO Recipe Execution Tool, version 3.7.2. The 
standard recipes were used to optimally extract each spec- 
trum. A ThAr arc lamp spectrum was used to wavelength 
calibrate the spectra, flux calibration was achieved using 
observations of standard stars. We then barycentrically cor- 
rected the wavelength scales of each of the spectra. The 2001 
and 2002 spectra were obtained using UVES in a dichroic 
mode as part of the ES O Supernova la Progenitor Survey 
|Napiwotzki et aLll200ll ). They cover a wavelength range of 
~ 3500 to 6650A with ~ 80A gaps at 4580A and 5640A, 
the spectral resolution is around R ~ 18500. The 2007 spec- 
tra cover a different wavelength range of 3300-4520A in the 
blue arm and 5680-7515A and 7660-9440A in the red arm 
thus Ha is not covered. These spectra have a resolution of 
R ~ 70000. 

We determined the orbital p hases of the 2001 a n d 200 2 
spectra using the ep hemeris of lO'Donoghue et all l|2003h . 
iParsons et al.l (|201dl ) showed that the eclipse times of QS 
Vir vary substantially after 2002 hence for the 2007 spectra 
we used the following linear ephemeris fitted to the eclipse 
times after 2006: 

MJD(BTDB) = 48689.135 529(54) + 0.150 757 6116(14)^. 



A Stellar Prominence in QS Vir 3 




4100 
Wavelength (A) 



5100 5200 
Wavelength (A) 




6100 6200 
Wavelength (A) 



Figure 1. 2001 UVES spectrum of QS Vir with the prominence features labelled. The sharp dip in the bottom panel labelled "inst" is 
an instrumental feature. Hydrogen Balmer, Ca II and Si I emission from the red star is also visible (the silicon feature is labelled) as well 
as many absorption features from the red star. There is also weak Mg II 448lA absorption though this is not visible on this scale. 



absorption is material within a stellar prominence originat- 
ing at the red star. Stellar prominences have been seen 
in a number of different systems; slingshot prominences 
have been observed in the cataclysmic vari ables BV Cen 
d Watson et all 120071 ). IP Peg and SS Cyg (|Steeghs et all 
1996) and long-lived loop prominen ces have been detecte d 
on the seconda r y star in AM Her (|Gansicke et al.l Il99ct ) . 
iDunstone et al.l (|2006l ) found a prominence from the K3 
dwarf Speedy Mic ( BO Mic) which shows similar features to 
those seen here and lPeterson et all (|20 ldT ) discovered a large 
coronal loop in the Algol system. Since the red st ar in QS 
Vir is known to be active (|Q'Donoghue et al1l2003l ). promi- 
nences originating at this star are certainly possible. These 
are the first observations of prominence material in absorp- 
tion within a binary; all of the previous examples show sharp 
emission features from prominence material within binaries. 
This provides useful diagnostics of the physical conditions 
within the prominence. 



4 MODELLING OF THE PROMINENCE 

The ratio of the equivalent widths (EWs) of the various 
absorption lines arising from the prominence can be used to 
constrain its temperature and density. S ince the red star in 
QS V ir is close to filling its Roche lobe (jO'Donoghue et all 
2003) these kind of features are important for understanding 
any accretion processes taking place in the system. 

As we remarked earlier the depth of the prominence 
absorption in the higher Balmer series is far deeper than 
the total light emitted by the red star at these shorter 
wavelengths hence the prominence material must be block- 
ing light from the white dwarf. However, in order to de- 
termine accurate EWs, the contribution from the red star 
has to be removed. We determined the red star's contri- 
bution to the flux across the entire wavelength range cov- 
ered by using a model spectrum of a 14,000K logg = 8.25 
white dwarf (calculat ed using TLUSTY and SYNSPEC; 
iHubenv fc Land Il995l) and a range of M dwarf template 
spectra ( Rebassa-Mansergas et alj|2007l ) from spectral type 
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Figure 2. The 2001 spectrum with the red star's contribution subtracted. The grey line is the M3 template that was subtracted from 
the original spectrum. This model starts at a wavelength of 3830A hence before this point we extrapolate the continuum level, since the 
white dwarf dominates at these wavelengths this has a minimal effect. The red star features have not been completely removed since the 
resolution of the template is less than the UVES spectra resolution but the overall flux level from the red star has been removed. The 
emission features were not included in the model and therefore have not been removed either. The small contribution from the white 
dwarf around the Ho absorption makes determining an accurate equivalent width difficult. 



M3 to M4 in order to clos ely match the parameters found by 
lO'Donoghue et al.l <|2003h . We used these to derive synthetic 
fluxes for both of the stars and adjusted th e overall level 
to ma tch the relative magnitudes shown in IParsons et alj 
(|2010h for the ULTRACAM u , g', r' and i' filters. We fit- 
ted the ratio of the in-eclipse and out-of-eclipse flux levels 
for the g' and r' bands since these best cover the Balmer 
series. There was little difference between the spectral types 
used for the secondary, however, the M3 model combined 
with the white dwarf model best fitted all of the filter ra- 
tios. Figure [5] shows the 2001 spectra with the M3 model, 
corrected to the radial velocity of the red star, subtracted 
off. Also shown is the M3 model which gives an indication of 
the relative contributions of the two stars at different wave- 
lengths. 

For the H<5 to Ha lines, the absorption from the promi- 
nence occurs on top of the emission from the secondary star 



because the emission becomes wider in these lines. Addi- 
tionally, the He absorption lies on top of the Ca II emis- 
sion from the red star. Ideally, in order to correct the EWs 
for this emission, we would subtract off the emission using 
observations of the lines at a similar phase without the ab- 
sorption component. However, since the UVES spectra cover 
only three different orbital phases, none of which are close to 
each other, we cannot subtract the emission from the spectra 
accurately, and therefore we use the white dwarf model to 
determine the continuum level and correct the EWs accord- 
ingly. This approach works well for the He to H/3 lines where 
the absorption occurs in the wings of the emission, but un- 
fortunately the Ha absorption lies right on the strong emis- 
sion peak and an accurate correction would require knowl- 
edge of the emission profile. Since this is unknown, we cor- 
rect the absorption to the white dwarf continuum but note 
that there is likely to be an additional correction required 
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Table 2. Measured parameters of the prominence absorption 
lines. "C Fac" is the additional equivalent width needed to correct 
for the emission from the secondary star. 



Line ID Velocity FWHM EW C Fac 

kms -1 kms -1 (mA) (mA) 

Ha 6562.760 -114.1 ± 1.3 74.0 ± 4.2 389 ± 30 600 

H/3 4861.327 -119.4 ±0.6 40.5 ± 1.5 604 ± 47 160 

H7 4340.465 -125.0 ±0.4 35.5 ± 1.1 506 ± 28 60 

H5 4101.735 -120.9 ±0.3 31.9 ±0.9 455 ± 32 15 

He 3970.074 -126.9 ±0.4 38.2 ± 1.3 400 ± 29 5 

Can 3968.469 -120.0 ±2.0 28.6 ± 5.3 89 ± 26 

Can 3933.663 -121.7 ±0.9 18.0 ± 2.3 82 ± 20 

H8 3889.055 -128.6 ±0.6 30.1 ± 1.5 322 ± 28 

H9 3835.397 -121.4 ±0.7 28.6 ± 1.8 232 ± 25 

H10 3797.910 -123.0 ±0.8 25.9 ±1.8 167 ± 25 

Hll 3770.634 -121.0 ± 1.0 28.2 ± 2.4 159 ± 22 

H12 3750.152 -130.8 ±1.8 29.3 ± 4.4 125 ± 23 

H13 3734.372 -124.9 ± 1.4 25.7 ±3.5 89 ± 19 

H14 3721.948 -127.9 ± 1.5 28.1 ±3.6 98 ± 17 

H15 3711.971 -118.6 ±2.6 33.8 ± 6.5 86 ± 19 



hence all measurements of this line will be unreliable. Addi- 
tionally, as Figure [2] shows, the contribution from the white 
dwarf around Ha is much lower than the red star's contri- 
bution making the absorption feature weaker still. 

4.1 EWs of the absorption lines 

We measured the EW, velocity offset and full-width-half- 
maximum (FWHM) of all of the prominence features. We 
determined the velocity offset and FWHM by fitting the 
line with a combination of a Gaussian and a straight line. 
The results of this are listed in Tabled 

After correcting for the red star emission, we find that, 
within the uncertainties, the fluxes in the centres of the 
sharp absorptions of H5, H7 and H/3 are zero thus the promi- 
nence must be large enough to block our view of the entire 
white dwarf, and also the material in the prominence cannot 
be very hot. A weighted average of the velocities of the lines 
gives the radial velocity of the absorption as —123.3 ± 0.2 
kms -1 . The radial ve locity of the white dwarf for this phase 
is — 118 ± 11 kms -1 (|Q'Donoghue et ai1l2003l . systemic ve- 
locity included). As previously mentioned, this is evidence of 
the absorption being from a prominence since any material 
in solid body rotation within the system, along the line of 
sight of the white dwarf, must have the same radial velocity 
as the white dwarf. 

4.2 Column densities 

In this section we describe the model we used to determine 
the column densities of hydrogen and calcium in the promi- 
nence and constrain its temperature and turbulent velocity. 

4.2.1 The model 

We calculate the line absorption profile using a Voigt profile 
of the form: 

m e c 2 yA\ D ^/TT 7rAA 2 +7 2 y 



where e is the electron charge, Ao is the central wavelength 
of the line, / is the oscillator strength, m e is the mass of 
an electron, c is the speed of light, the star symbol is the 
convolution operator and A A = A — Ao, the distance from 
the centre of the line. The first term inside the brackets is a 
Gaussian profile representing the thermal and microturbu- 
lent broadening, in this case 

AA D ^(^-Kf 5 , (3) 
c \ m J 

where k is the Boltzmann constant, T is the temperature, 
m is the mass of the atomic species involved and £ is the 
microturbulence. This is convolved with the Lorentzian pro- 
file on the right in Eqn[2]which represents the natural line 
broadening where 7 is the FWHM given by 

where the Aij are the Einstein A-coefficients contributing to 
the line. This gives us the absorption profile per atom. 

The optical depth is then given by multiplying this by 
the column density (N - the total column density of the 
atomic species in a given energy level): 

T\ = ot\N. 

Hence integrating over the line profile gives us the EW via 

EW = j(l- e~ Tx )d\. (5) 

Therefore the only variables are the temperature, T, the 
turbulent velocity, £, and the column density, N. However, 
there is a degeneracy between the thermal and turbulent 
broadening whereby Vx urb + Vx hcrm = C, a constant. For 
large background sources the EW is also dependent upon the 
covering fraction (the amount of the emission source covered 
by the absorber), which leads to a degeneracy between the 
column density and the covering fraction. Fortunately, since 
the white dwarf is small it is completely covered (as can also 
be seen in the very deep absorption cores in Figure [2]| hence 
the degeneracy is broken and we are able to obtain reliable 
column densities. 



4-2.2 Column density of hydrogen in the n = 2 level 

We fitted our model to the measured EWs of the hydro- 
gen Balmer absorption lines (with the ex ception of Ha) us - 
ing Levenberg-Marquardt minimisation (|Press et al.l 1198(f ) 
for a range of turbulent velocities. We find that N2 = 
1.26 ± 0.10 x 10 14 cm -2 . For temperatures below 7600K the 
microturbulence can balance the thermal broadening thus 
offering no constraint. The fit to the EWs for a temperature 
of 7600K is shown in Figure [3] By using a range of turbulent 
velocities we also find that V^ ulh + VT hcrm = 63.2 (kms -1 ) 2 . 
The three sigma upper limit on the temperature from fit- 
ting the EWs is 9950K, as also shown in Figure [3] Above 
this temperature the lower lines of the series become unsat- 
urated and the EWs are overestimated. However, these high 
temperatures (above 7600K) require zero turbulence to fit. 
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Figure 3. Model fits to the EWs of the Balmer series. 'V is the 
upper energy level involved in the transition (Ha = 3, H/3 = 4 
etc). For temperatures below 7600K the EWs can be fitted by in- 
creasing the turbulence. The 9950K model shows the three sigma 
upper limit to fitting the EWs, any additional increase in temper- 
ature causes the lower lines of the series to become unsaturated. 
Hence a good fit is still obtained if there is no turbulence and a 
high temperature (as is favoured by ionisation equilibrium con- 
siderations). Due to the uncertainty in the correction required to 
the Ha line due to emission from the red star we do not fit this 
line. The dashed line is J r\d\ which shows the predicted EW 
without the effects of saturation; the higher lines of the series 
show very little saturation. 



4-2.3 Column density of Ca II 

Since we also detect the calcium H and K lines from the 
prominence we can carry out the same analysis as for the 
hydrogen Balmer lines in order to establish the column den- 
sity of Ca II. Unfortunately, since we have only two lines with 
rather large uncertainties in the EWs we cannot constrain 
the column density as tightly as with the Balmer lines. 

We calculated the range of column densities that will 
fit the observed EWs of the Ca II lines (to within 2a) 
over a grid of temperatures from 3000-1 1,000K and tur- 
bulent velocities from 0-8 kms -1 . We find that high tur- 
bulent velocities (> 4 kms -1 ) favour low column densities 
(11.3 < log 10 -/Vcaii(cm -2 ) < 11.6) for all temperatures. 
However, lower turbulence leads to a much larger spread in 
column densities (11.5 < log 10 A r caii(cm^ 2 ) < 14.6) with a 
much larger dependence on temperature whereby the lower 
the temperature, the higher the column density. 

4-2.4 Column density of hydrogen in the n = 1 level 

We can constrain the column density of neutral hydrogen 
using our constraints on the column densities of hydro- 
gen in the n — 2 state and calcium. We assume a size of 
the prominence of 10 10 cm (see section l5~Tj) and solar abun- 
dances. From the Saha-Boltzmann equation (assuming local 
thermodynamical equilibrium) we can determine the column 
density of neutral hydrogen that would give the observed col- 
umn density of calcium. We can also use the column density 
of hydrogen in the n = 2 state to calculate the column den- 
sity of neutral hydrogen via a Boltzmann factor. The results 
of both of these calculations are shown in Figure [4] over a 
range of temperatures. The large uncertainty in the column 
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Figure 4. Column density of hydrogen in the prominence as a 
function of temperature. Using the Ca II lines permits a large 
range of column densities over the temperature range. Using the 
Balmer lines places tighter constraints on the column density. The 
dashed line shows the three sigma upper temperature constraint 
from fitting the Balmer lines. The permitted column density of 
hydrogen from the Balmer lines and Ca II lines overlap at a higher 
temperature (~ 8800K) which is still permitted from fitting the 
Balmer lines but requires very low turbulence. 



density of calcium results in a large range of possible neu- 
tral hydrogen column densities. However, the hydrogen in 
the n = 2 state gives a much tighter constraint. The two 
regions overlap at temperatures above 8800K implying that 
the prominence material is hot. This temperature is below 
the three sigma limit set by the EWs of the Balmer lines 
and implies low levels of turbulence (as seen in Figure |3J. It 
is also possible that the material may not be entirely in lo- 
cal thermodynamical equilibrium. However, our various con- 
straints are consistent with the prominence having little or 
no turbulence and a temperature of ~ 9000K. These results 
are consistent with obs ervations of other stellar prominences 
Dunstone et alj 20061) and modelling of solar prom inences 



Gilbert et al.ll2006l ; lLabrosse fc Gouttebrozdl200ll ). 



5 DISCUSSION 

5.1 The size of the prominence 

We can place a lower limit on the size of the prominence 
since we know that the white dwarf is completely covered 
hence the prominence must be at least as large as the white 
dwarf. However, we can do better than this since we know 
that the white dwarf is obscured for the full 300 second expo- 
sure therefore the minimum size of the prominence must be 
2nA4>R — 0.145-R, where A(f> is the change in phase during 
the exposure and R is the distance of the prominence from 
the white dwarf. The width of the narrowest feature tells us 
that R > 0.22a (see section [372|) hence the size of the promi- 
nence must be at least 0.22a x 0.145 = 0.032a ~ 3.7i?wD- It 
could of course be considerably larger than this limit. 

A minimum distance of the prominence from the red 
star can be calculated since we know that the prominence 
lies along the line of sight from the white dwarf at an 
orbital phase of 0.16. This gives a minimum distance of 
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Figure 5. Left: a visualisation of the system as viewed from above. The dotted lines show the Roche lobes of the two stars and the 
path that a stream would take were one present. Also marked are the positions of the L4 and L5 points, the two straight lines show the 
viewing angle of the white dwarf from Earth at the start and end of the exposure, thus the prominence material is located within the 
shaded a rea. Right: a velocity ma p of QS Vir with the same features as the left panel. The star shows the location of the emission feature 
found bv lO'Donoghue et al. the circle represents the uncertainty in its position. 



the prominence fr o m th e centre of the red star of 0.81a. 
lO'Donoehue et al.l (|2003l ) estimate i? RD = O.42i? which 
gives a minimum distance of the prominence above the sur- 
face of the red star as 0.62i?Q, almost 1.5 stellar radii away. 

Limits can be placed on the density of the material 
within the prominence. From section 14.2.41 we determined 
a range of column densities of neutral hydrogen of 0.5 x 
10 22 < Nh (cm -2 ) < 1.0 x 10 22 . Using our minimum promi- 
nence size (0.032a) gives an upper limit on the density of 
nu = 3.5 x 10 12 cm -3 . A loose lower limit can be calculated 
by assuming that the prominence is smaller than the binary 
separation, this gives a density of nu = 5.6 x 10 10 cm~ 3 . 
This range of densities is consistent with observations and 
modelling of both solar Jvan Ballegooiien fc Cranmerll201ol ) 
and other stellar prominences (|Schroederlll983l ). 

The fact that QS Vir is an eclipsing system and we 
know that the prominence passes in front of the white dwarf 
means that it must also pass in front of the red star. There 
is no evidence in the 2002 or 2007 spectra of any prominence 
features; this could either be due to the prominence being 
a short-lived feature or that we have not covered the phase 
at which the prominence crosses the face of the red star. 
Observations at a range of phases including the detection of 
the orbital phase at which the prominence passes in front of 
the red star would tightly constrain the size, location and 
mass of the prominence. 



5.2 Evidence for a detached system 

Figure [S] shows a visualisation of the system. All the features 
were drawn using the bes t estimate system parameters from 
lO'Donoehue et alj (|2003l ). we used a mass ratio of q = 0.52 
and K WD + K Bec = 137 + 266 = 403 kms" 1 . We note here 
that this is not an exact visualisation, due to the uncertainty 
in the system parameters, but it serves to understand better 



the location of the prominence and other features previously 
identified in this system. 

The left hand panel of Figure [5] shows the system in 
spatial co-ordinates, we also show the positions of the Roche 
lobes of the two stars and the route that a stream would 
take (were one present). The shaded region indicates the 
potential positions of the prominence (it must lie along the 
line of sight of the white dwarf at <f> = 0.16). Interestingly, it 
is possible for the prominence to be located close to the L5 
point. Since this is an unstable, though neutral point, the 
material would still have to be supported in some way (e.g. 
by a magneti c field) . Similar features have been identified by 
IJensen et all (|l986t ) in the eclipsing PCEB V471 Tau. They 
found dips in the soft X-ray flux from material along the line 
of sight at orbital phases consistent with the material being 
at both the L4 and L5 points. They attribute this material to 
very l arge coronal loops anchored at the K star. [Wheatle y] 
( 1998) also found X-ray dips in V471 Tau but caused by 
material much closer to the K star (within a stellar radius). 

The right hand panel of Figure[S]shows the system in ve- 
locity co-ordinates. The same features as the left hand panel 
are shown (Roche lobes and stream) as well as the possible 
positions of the prominence (assuming that it is locked to 
the red star and hence the binary). We also indicate the 
p osition of an emissi o n fea ture seen in the Doppler image 
of lO'Donoghue et all (|2003l ) at (V x ,V y ) = (75, -75) kms" 1 
the circle around it is the uncertainty in its position given 
the resolutio n of their Doppler i mage. This feature was in- 
terpreted bv lO'Donoghue et al.l l|2003l ) as emission from an 
accretion stream, however a stream would not be located 
at these co-ordinates in velocity space, rather it would be 
located outside the Roche lobe of the white dwarf along 
the nearly diagonal line to the upper-right of the primary's 
Roche lobe as sh own in Figure ^ Therefore the emission fea- 
ture detected bv lO'Donoghue et alj <|2003l ) cannot be from 
an accretion stream. 
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The range of possible positions for the prominence 
passes through the emission feature, so, rather than being 
due to a low level accretion stream shock, the emission fea- 
ture seen in the Doppler image could well be related to a 
prominence. If the prominence we have detected is due to the 
same material as seen in the Doppler map then it is remark- 
ably lo ng lived since the spectra used to create the Doppler 
map in lO'Donoghue et all l|2003h were taken between April 
1992 and June 1993 almost a decade before the UVES spec- 
tra presented here, however, further data are needed to test 
this. 

These observat ions support the argument made by 
[Rib eiro et all (|2010l ) that the material physically inside the 
white dwarf's Roche lobe is likely to be the result of wind 
accretion onto the white dwarf. Their claim that this may 
be facilitated by a prominence-like magnetic loop appears to 
be supp orted by our data. O ur observations, combined with 
those of lRibeiro et al.l (|2010l ). show that QS Vir is currently 
a detached system. 

Mg n 448lA absorption is present in all of the spectra, 
though it is quite weak in the 2001 and 2002 spectra. The 
strength of the absorption in all three spectra are consistent 
within their uncertainties. The feature is best seen in the 
2007 data due to the longer exposure time and is shown in 
Figure [5] Since this absorption is seen in all the spectra it 
cannot have the same origin as the sharp features seen only 
in the 2001 spectra. Furthermore, the velocity of the line 
is in the opposite sense to the red star meaning that this 
absorption comes from the atmosphere of the white dwarf. 
This feature potentially allows us to measure the radial ve- 
locity of the white dwarf, however, with just three spectra a 
measurement of the radial velocity is not possible with the 
data to hand. QS Vir's current evolutionary stage (pre-CV 
or hibernating CV) can also be addressed by our spectra; 
the identification of Mg II 448lA absorption from the white 
dwarf allows us to constrain the rotational velocity of the 
white dwarf. 

lO'Donoghue et al.l (|2003h determined a rotational ve- 
locity for the white dwarf of V TO t = 400 ± 100 kms" 1 from 
modelling of Hubble Space Teles cope /STIS (i/ST/STIS) UV 
spectra, suggesting that the white dwarf had been spun-up 
by an earlier phase of mass transfer. However, this measure- 
ment is marginal since it is at the limit of the HST/STIS 
spectral resolution. Figure [S] shows the Mg II 448lA absorp- 
tion feature in the 2007 spectrum with three models with 
different v sin i values fo r the white dwarf over plotted (cal- 
culated using TLUSTY; iHubenv fc Lan3ll995l ). The model 
with usinz = 400 kms -1 overestimates the width of the 
line and even with a rotational velocity of 100 kms" 1 the 
predicted line width is too wide, a rotational velocity of 
vsini = 55 kms" 1 fits the line best. However, the white 
dwarf moves by ~ 100 kms - during the 1800 second expo- 
sure which can entirely explain the width of the line, mean- 
ing that a v sin i = is favoured. Nevertheless, we can rule 
out a lar ge vsini value, wh i ch wa s one of the main argu- 
ments of lO'Donoghue et al.1 (|2003l ) that QS Vir is a hiber- 
nating CV. Hence, the fact that there is no evidence of mass 
transfer in the system, leads us to conclude that QS Vir has 
yet to reach a semi-detached configuration and is therefore 
a pre-CV rather than a hibernating CV. 
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Figure 6. Mg II 448lA absorption from the white dwarf seen 
in 2007 in slightly poor conditions. Overplotted are three mod- 
els with varying vsini values of 55 kms -1 , 100 kms -1 and 400 
kms -1 . The vsini = 55 kms -1 model fits the line best but the 
quality of the data is not good enough to tightly constrain its 
value. However, the narrow nature of this feature confirms that 
the white dwarf is not a rapid rotator. 



6 CONCLUSIONS 

The main sequence star in the eclip sing PCEB QS Vir is 
very c lose to filling its Roche lobe. ISchreiber fc Gansickel 
( 2003) showed that QS Vir is one of the few PCEBs that has 
passed through most of its PCEB life time, this makes it a 
particularly interesting system since the majority of PCEBs 
are young systems that have only just emerged from their 

com mon envelope phas e . 

lO'Donoghue et all (|2003l ) found an emission feature 
from within the Roche lobe of the white dwarf which, com- 
bined with their measurements of a rapidly rotating white 
dwarf, led them to believe that QS Vir was in fact a hi- 
bernating CV that had become detached as a result of a 
nova explosion. Here we have used UVES spectroscopy to 
identify and characterise a stellar prominence in the sys- 
tem w hich lies along the line of s ight of the emission feature 
seen in lO'Donoghue et alj (|2003l ) leading us to conclude that 
the emission feature may be associated with a prominence 
rather than a low-level accretion stream. It is also possible 
that the prominence may be located at the L5 point of the 
system. We have also detected Mg II absorption arising from 
the atmosphere of the white dwarf; the width of this feature 
is incompatible with rapid rotation of the white dwarf sug- 
gesting that QS Vir has yet to become a CV and is thus a 
pre-CV. 
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